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Abstract. We study the full probability distribution of the charge transmitted through a mesoscopic diffu- 
sive conductor during a measurement time A. We have considered a semi-classical model, with an exclusion 
principle in a discretized single-particle phase-space. In the large A limit, numerical simulations show a 
universal probability distribution which agrees very well with the quantum mechanical prediction of Lee, 
Levitov and Yakovets [PRB 51 4079 (1995)] for the charge counting statistics. Special attention is given 
to its third cumulant, including an analysis of finite size effects and of some experimental constraints for 
its accurate measurement. 

PACS. 72.70.4-m Electronic transport in condensed matter: Noise processes and phenomena: Electronic 
transport in mesoscopic or nanoscale materials and structures - 73.50.Td Electronic structure and electrical 
properties of surfaces, interfaces, thin films, and low-dimensional structures:Noise processes and phenomena 
- 71.10.Fd Electronic structure of bulk materials:Lattice fermion models 
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1 Introduction 

The electrical current across a conductor results from a 
flow of charge carriers. Their discreetness is responsible 
for current fluctuations Si -call ed shot no ise- character- 
ized by a noise power Si = 2 J Si(t)Si(0)dt proportional 
to the mean current I. For example, a current I of uncorre- 
cted charge e will exhibit a full shot noise Si = 2.e.I. For 
a recent review on shot noise and for general references, 
see|. 

In diffusive conductors shorter that the electron-phonon 
mean free length, an universal shot noise power Si = 
(2.eJ)/3 has been predicte d and vali- 

dated experimentally 0,0[L2]|]^. 

This 1/3 reduction over the full shot noise focused a 
lot of attention because it has been derived in two differ- 
ent frameworks: quantum mechanics (scattering matrix 
theory and Green function technique fj,^]) and 

semi-classical mechanics (Boltzmann-Langevin approach 
including Pauli's exclusion principle @,|§|,f§)- While Lan- 
dauer saw a numerical coincidence between these results Jl4j 
de Jong and Beenakker minimized the role of phase 
coherence, which entails that an exclusion principle is the 
only irreducible concept behind the two treatments. 

To be more complete, we should also mention the stud- 
ies of the shot noise reduction in non-degenerate diffusive 
conductors: @,0,|8|. 

The shot noise power can be seen as the statistical 
variance of the measurable current fluctuations, (by mea- 
surable current, we mean: current averaged over a time 
scale much larger than the electron-electron correlation 



time, which is always the case in electrical set-ups, even 
for very highpass systems.) In addition to the variance, the 
shot noise should also be characterized by all the higher 
moments of current fluctuations. In this paper, we will fo- 
cus on the information beyond shot noise power, in the 
perspective of testing the merits of classical versus quan- 
tum mechanics in the noise reduction process. 

The first significant result beyond the shot noise power 
was derived by Levitov and Lesovik [|l9| . Starting from the 
description of a charge-counting operator, they derived the 
full statistics of transmitted charges through a point con- 
tact, or in other words the probability that Q charges are 
transmitted through a lumped scatterer during a given 
measurement duration A that is large compared to the 
electronic correlation times. Note that in this large A 
limit, Q/A become what we called a measurable current. 



Two years later, Lee et al. [20 combined this result with 



the bimodal distribution of the transmission eigenvalues 
predicted by the random matrix theory for a phase co- 
herent diffusive conductor |2lj. In the large A limit, they 
obtained the statistics P q (Q) for the transmitted charges 
Q by this fully quantum derivation. 

In this paper, we will show that the full statistics of 
transmitted charge through a diffusive conductor is a semi- 
classical quantity. To do so, we show that the prediction of 
Lee et al. can be recovered in a semi-classical model which 
accounts for an exclusion principle but not for phase co- 
herence. This result generalizes to diffusive conductors an 
earlier derivation for double-barrier junctions j23|. 
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2 The model 

A minimal modeling of out-of-equilibrium semi-classical 
mesoscopic conductor consists in a ID open system with 
a flow and a back-flow of charges of opposite velocities. 
This degenerate model can be easily adapted to a one 
dimensional chain: each site being associated with a back- 
scattering probability and a Pauli exclusion rule. Such de- 
generate systems are often considered in the theoretical 
physics literature as "simple exclusion process" or traffic 
models. Indeed some of the analytical results on shot noise 
power mentioned previously have been derived with such 
models |8[p^] . In order to strictly overlap with these pre- 
vious studies, we stick with the exact sequence of dynamic 
rules considered by Liu, Eastman and Yamamoto |2^j . 

The model consists in a chain of N sites. Each site is ei- 
ther empty, occupied by a right or a left-moving charge or 
by two charges propagating in opposite directions. During 
each time step, all charges are shifted to the next site, ac- 
cording to their direction. After this, each charge is likely 
to back-scatter with a probability r, provided that the 
resulting state is empty. The system is maintained out of 
equilibrium by dissymmetrical boundary conditions: charges 
are injected at each time step at one end of the chain while 
both ends act as perfect absorber for charges incoming 
from the chain. 

For the model parameters considered in this paper, 
over 10 s samples are sometimes required for the statistics 
to converge up to the desired accuracy. In practice, the nu- 
merical simulation could still be performed on a desktop 
computer thanks to a code which core is restricted to low 
level-processor-instructions. Basically, each configuration 
of the chain occupation can be coded by the binary rep- 
resentation of two integers, one for each charge direction. 
Each time step is a combination of a register shift (right 
or left), a bit increment to account for the charge injec- 
tion at the boundary, and a bit-to-bit comparison of the 
two integers to check for the scattering which are compat- 
ible with the Pauli exclusion principle. The effectiveness 
of scattering is set by precomputed series of random bits, 
refreshed when necessary. A special attention was dedi- 
cated to the validation of the random numbers algorithm. 
In addition to this direct simulation, two semi-analytical 
methods provided a cross validation of the results up to 
N=6 sites. For longer chains, those methods were either 
too CPU-time or memory consuming. 



3 Results 

The dots on Figure 1 present the probability distribu- 
tion of the transmitted charges Qa for typical parame- 
ters of the semi-classical model (measurement duration: 
A = 18000 time units [p4| , site's back-scattering prob- 
ability: r = 50%). For comparison, the distribution Q q 
predicted by Lee et al. quantum model is plotted for the 
same average transmission (thin line) pp| . 

The central limit theorem states that both these distri- 
butions converge to a gaussian distribution when A — > oo. 
In this large A limit and within finite size corrections, it 
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Fig. 1. Left scale: probability distribution from the number of 
transmitted charges during A = 18000 time units. Each site 
presents a r — 50% back-scattering probability, which results 
in a G = (Qa)/A — 9% conductivity. Dots: present semi- 
classical model. Thin line: Lee et al. quantum model for the 
same average transmission E3]. Right scale: difference between 
each distribution and their gaussian fit. Open circles: present 
model, thick line: Lee et al. model. 



has been demonstrated that the distributions variances of 
both models -i.e. the shot noise power- become equal 
fjOf . More precisely, in the large A limit, de Jong and 
Beenakker showed that when the number of sites N in- 
creases, the Fano factor ((Qa — (Qa)) 2 )/(Qa) converges 
to ((Q q ~ (Q q )) 2 )/(Q q ) = 1/3. _ 

Beyond the gaussian approximation, the non-gaussian 
contribution of both distributions is also plotted on fig- 
ure 1 (right scale). The central finding of this work is the 
excellent agreement between the semi-classical (open cir- 
cles) and the quantum models (thick line) which strongly 
suggests that a semi-classical picture fully accounts for 
the whole statistics of transmitted charges in a diffusive 
conductor. 



4 Finite size effects: the signature of 
correlations 

The remaining of this paper will focus on the finite size 
effects associated with the three parameters of the model: 
the chain length or number of sites N, the measurement 
duration A and the back-scattering probability r on each 
site. 

The third cumulant ((Qa — (Qa)) 3 ) provides a useful 
measure of the deviation from the gaussian distribution 
and it can be directly compared with the quantum model 
prediction ((Q q - (Q q )) 3 ) = 1/15. (Q q ) |3. In this last re- 
lation, the linear scaling with (Q q ) is simply a consequence 
of the fact that for large measurement duration (A = oo 
in this model), (Q q ) is the sum of almost uncorrelated 
random variables which effective number scales linearly 
with A. One of the motivation for considering cumulants 
is precisely that they behave linearly with respect to the 
addition of independent variables. Therefore, the physics 
of the short time scales electron correlations is captured 
by the numerical prefactor 1/15. In the following, it will 
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Fig. 2. Dependence of the third Fano factor F3 = {(Qa — 
(Qa)) 3 )/(Qa} with the measurement duration A for the 
present model (solid lines) and for the same model without 
exclusion principle (dashed line). For all curves the back- 
scattering probability is r = 85%. The number of sites per 
chain is <0: 2, triangles: 4, circles: 7. The transient regime evi- 
dences an anti-correlation time between the charge arrivals. 



be useful to call the ratio of the third cumulant by the av- 
erage transmission the third Fano factor F3, in reference 
to the usual Fano factor F. 

Figure 2 presents the semiclassical third Fano factor 
Fa = ((Qa — (Qa)) 3 }/{Qa) versus the measurement dura- 
tion A for different combinations of N and r (solid lines) . 
After a transient regime, the signal settles to a constant 
level. The transient regime duration defines a correlation 
time r. Beyond this regime, for A >> t, the random vari- 
able Qa scales like A/t as expected and the third Fano 
factor tends towards a constant value. Note that the mea- 
surement duration A — 18000 used in figure 1 -to compare 
Qa and Q q - is validated by a plot similar to figure 2 for 
r = 50%. 

The simulations show that r increases when the chain 
length is increased as one would expect for a similar model 
without exclusion principle. In this latter case, we shall 
show that t is of the order of the scattering time through 
the chain. Let us compare further the model with and 
without exclusion principle. On figure 2, F3 versus A is 
plotted for this second model for the same back-scattering 
probability and chain lengths (dashed lines). The transient 
regimes are significantly larger in this second model. In the 
A = limit, F3 = 1 as expected from a Poisson distribu- 
tion while for large A, F3 converges to (1 — G).(l — 2G) 
where G = (Qa)/A is the conductivity. More extensive 
calculations confirmed this large A limit, which corre- 
spond to a binomial distribution. 

This result can be understood with the basic sketch of 
figure 3. To begin with, the main feature of the model in 
the absence of exclusion principle is that electrons emitted 
at different times are completely uncorrelated. The prob- 
ability distribution of the transmitted charge is therefore 
completely determined by the probability distribution for 
a single particle to leave the system on the right, as a 
function of the time spent inside it. Let us denote by r 
the width of this distribution. The left picture on figure 
3 illustrates what happens as A 3> r. Most detected par- 



ticles have been emitted in region 2. For them, the mea- 
surement window covers the whole span of likely values for 
their lifetime in the system, before leaving on the right. 
Region 1 corresponds to particles emitted early, so they 
have to spend more time than on average in the system 
in order to contribute to the signal. Symmetrically, region 
3 corresponds to particles spending a shorter time than 
on average in the system. Region 1 and 3 have a width of 
order t, whereas region 2 lasts for a time equivalent to A 
as A 3> t. All the particles emitted in region 2 contribute 
independently to the signal with the same probability , 
so the measured charge follows a binomial distribution: 
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The right picture on figure 3 illustrates the opposite 
situation, when A <C r. First, we expect anticorrelations 
between charges measured in the time intervals [t, t + A] 
and [t',t' + A] if A <| t - t' |< r, simply because both 
measurements involve particles which have been emitted 
in a same time interval of width proportional to r, and 
the same particle cannot be detected twice. As 4 < r, 
the probability for each particle to contribute to the sig- 
nal becomes very small, so we get a Poisson distribution. 
The cross-over from this Poisson distribution to the bi- 
nomial one occurs for A ~ r, and is precisely due to the 
correlations between the small time intervals [t, t + A] and 
[t',t' + A] when A <| * - t' |< r. The main value of 
this overly simple model is to illustrate the role of the 
correlation time r. It also emphasizes the need to give 
an accurate modelling of the injection process. The bi- 
nomial law for A t is strongly connected to the fact 
that the particles are injected periodically in the system. 
If the injection process is instead taken to follow a Poisson 
distribution, the scale r doesn't appear and the transmit- 
ted charge obeys a Poisson distribution, for all values of 
A. These two features (role of r and importance of the 
injection mechanism) are of course expected to be found 
in more complex and more realistic models, such as the 
model simulated here, although full analytical treatments 
are not easily available. Let us now return to the case with 
the exclusion principle. 

Figure 4 presents the Fano factor F and the third Fano 
P3 factor versus the conductance G — (Qa)/A for vari- 
ous numbers of sites N. This plot completes a similar one 
published in |2^] for the Fano factor. The continuous lines 
are the factors F = (1 - G) and P 3 = (1 - G).(l - 2G) 
for a point scatterer predicted by the coherent scattering 
formalism ||(| . The agreement with the data in the G — ► 1 
limit can be understood easily: for N.r < 1, the mean free 
path becomes longer than the chain and the whole sys- 
tem can be considered as a lumped (or point) scatterer. 
On contrary, in the low conductivity limit, each charge 
undergoes many collisions in the chain and a diffusive- 
like behavior is expected. It is known that as the chain 
length increases, the data for the Fano factor converge to 
a plateau at 1/3 [p||20|. Figure 4 shows that the same type 
of asymptotic plateau emerges for the third Fano factor 
data. A plateau at P3 = 1/15, plotted in figure 4, is com- 
patible with the data. 
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Fig. 3. A simple way to analyze the model in the absence 
of exclusion principle. Electrons emitted at different times are 
uncorrelated. They have a finite probability to contribute to 
the signal if their typical arrival time (lying in the tilted band 
of slope unity and width r) lies in the measurement window 
(the horizontal band of width A). Left large A limit (4>t). 
Right: small A limit (A -C t). The picture shows two different 
measurement windows separated by an interval shorter than r, 
which leads to strong anticorrelations between the two signals. 
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Fig. 4. Dependence of the Fano F and third Fano factors F3 
versus the conductivity G = (Qa)/A. The symbols are the 
data for N=l, 2, 3, 6, 10, 28 sites ; the dotted lines through 
those symbols are guideline. The thin lines are the dependences 
for a lumped scatterer predicted by the coherent scattering 
formalism: (1 — G) and (1 — G).(l — 2G). The two plateaux at 
F — 1/3 and F3 = 1/15 are predicted by Lee et al.'s quantum 
derivation for a diffusive conductor. 



Another point of view on these data supports the con- 
jecture of an asymptotic plateau at 1/15. Figure 5 presents 
the third Fano F3 factor versus the number of sites N for 
various back-scattering probabilities r. The residual dis- 
crepancy between the 1/15 limit and the N = 10 data 
appears as a finite size effect on N. The inset shows the 
difference between the third Fano factor and 1/15 as a 
function of the number of sites N for r — 80%. It is in- 
teresting to note that the convergence is compatible with 
a power law with a (—1) exponent, which may be related 
to long range correlations between charge carriers induced 
by the exclusion principle. 
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Fig. 5. Dependence of the third Fano F3 factors with the 
number of sites N. The plateau at F3 = 1/15 is predicted 
by Lee et al. Inset: F3 — 1/15 versus N for r — 80% up to 
N = 28. 

5 Perspective for experiments 

To the knowledge of the authors, experiments on the statis- 
tics of the transmitted charge have never been extended 
beyond the shot noise power. This may not be surpris- 
ing since such experiments are difficult for a fundamental 
reason: the physical phenomena revealed by higher cu- 
mulants are related to the correlation time scale r while 
the experimental probes (amperemetre, . . . ) always have 
a time constant A such that A 3> r. As we recalled ear- 
lier, in this limit the statistics of transmitted charges are 
very close to a gaussian and a large number n of statistical 
samples are required to discern the non-gaussian contri- 
butions, such as the third Fano factor F 3 . In the following, 
we show that under specific experimental conditions, these 
contributions remain measurable. 

We shall first consider the experimental analog of the 
cumulants ((Qa — (Qa)) 2 ) and ({Qa — (Qa}) 3 ) from which 
F and F3 are derived. A central point is the relation be- 
tween the number of transmitted charges Qa and the cur- 
rents i measured by the amperemetre ; our analysis differs 
from a previous one on this point p5| . The dynamical re- 
sponse of the amperemetre is characterized by a cut-off 
frequency fg or the corresponding time scale A = 1//b 
which can be seen as the duration of the charge counting 
from which each current output is inferred. Consequently, 
the relation between Qa and each current output i should 
then be: 

Qa ~ A.i - i/f B 

At this point, two comments are necessary. Firstly, fs 
should be understood as the effective frequency at the 
amperemetre output, which is set by the whole set-up, 
including filters, cables,. . . Secondly, in practice, the mea- 
surement chain is rarely a low pass but rather a bandpass, 
in order to filter out the low frequency noise or as a conse- 
quence of the typical specifications of RF elements. In the 
following, fs should be understood as the measurement 
bandwidth rather than the cut-off frequency. Note that 
we are implicitly assuming that the shot noise is white, 
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which is reasonable for the experimental conditions con- 
sidered below. 

Using the previous equation, the Fano factors F and 
i*3 can be written as : 

(T^3 = F 3 l(e.f B ) 2 

where the overline denotes the average over the band- 
width f B . The experimental uncertainty on F3 is mostly 
due to the uncertainty on (i — i) 3 , which results itself from 
undesirable current noises %n in the circuit (such as am- 
plifier noise and thermal noises) and from the finiteness 
of the number of statistical samples n. Assuming that i 
and ijy are independent and quasi-gaussian, the resulting 
signal to noise S/N of (i — i) 3 and F3 is indeed: 



S/N = (i - i) 3 /y/var((i + i N ) 3 ) 

where we used the relation between the variance var of 
the statistical estimate of the third cumulant of a gaussian 
and its second cumulant. To a reasonable approximation 

for our purpose, [i — i) 2 is the sum of the shot noise 2F.e.i 
and the thermal noise Ak.T/R across the conductor of re- 
sistance R: 



(i - i) 2 ~ (2F.e.i + Ak.T/R).f B /2 
If Sn is the noise power density of ijv, we also have: 

iff = Sn-Zb/^ 

The four last equations can be combined to give n as 
a function of the other parameters including S/N and i. 
If the amperemetre output is sampled at a rate f s (over- 
sampling being excluded) and if we set S/N = 10, the 
duration T of the experiment will be: 

_ n ^(2F.e7i + Ak.T/R + S N ) 3 
T = -f ~ 75. j 

fs f s .f B /.e\F 2 

For f B = f s _ = 1GHz, S N = 10- 24 A 2 /Hz, T = 4JC, 
R = 50.!? and i = 10 — 1000/iA, we find an experiment 
duration T < lmin, comfortable for an experimentalist. 

The typical RF frequencies that we find for f B has con- 
sequences on the procedure for processing the ampereme- 
tre output Q . A calibrated analogic device can certainly 
perform the i — ► i 3 function but even if a fully analog sig- 
nal processing is feasible p7] , we believe that a numerical 
acquisition is preferable here. Given the typical param- 
eters f s and T, real time storing of all the data is not 
possible with present technologies but a histogram of the 
measured current is enough for a post-processing of the 
cumulants and real time histograms can be acquired at 
frequencies up to several GHz. This signal processing pro- 
cedure would also allow to extract more information than 
just F3. In particular, figure 1 shows that the deviation 



from the gaussian has a shape which is nearly or- 
thogonal to the gaussian and not efficiently sensed by the 
i — » i 3 projection function of the third cumulant. Since 
this function is known theoretically, it could be used 
as a projection vector for the experimental histograms. It 
is clear that such a processing would significantly increases 
the S/N ratio. 

6 Concluding remarks 

The main result of this work is the perfect quantitative 
agreement between the classical model (with an exclusion 
principle in phase space) and the full quantum treatment 
of reference fl2G] for the complete probability distribution 
of transmitted - ^ charge through a ID diffusive conductor, 
during a measurement time A. Such an equivalence has 
already been discussed by many researchers in the last 
decade, either at the level of the second moment of this 
distribution (the noise power) or for the special 

case of a double-barrier system p3fT However, our impres- 
sion is that a first principle understanding of why this 
should be true is still lacking. Indeed, an advantage of the 
direct numerical simulation we have performed is to yield 
to the complete stationary out of equilibrium probabil- 
ity distribution on the configuration space of our system, 
which is the analogue of the Sinai-Rucllc-Bowen measure 
for continuous dynamical systems [p8[ . This information 
goes in principle much beyond the one contained in the 
Boltzmann stationary distribution or even the more re- 
fined Boltzmann-Langevin approach to the fluctuations in 
the single particle distribution around its stationary value. 
The latter formalism has proved very effective for the com- 
putation of noise reduction factors in various models 
y . But it relies on a key assumption on the fluctuations of 
the density current in single particle phase-space, namely 
it is a Poissonian random process with local correlations 
in space and time. Although this is a very reasonable as- 
sumption, as argued in |29[| , a rigorous derivation of this 
behavior requires in principle the exact knowledge of two 
particle correlations from the stationary distribution in 
configuration space. This goes far beyond the knowledge 
of the Boltzmann distribution. A natural question arising 
from the present work is to compare these two particle 
correlations obtained in the numerical simulation with the 
assumption made in [^] . We hope to get a clear answer on 
this important point in the near future. In fact, the strong 
similarity between our model and the simple exclusion 
process considered in |}0| suggests that the stationary dis- 
tribution in configuration space is a very complex object. 
Recently, Derrida et al have uncovered some fascinating 
properties of the rare large fluctuations of such distribu- 
tions, which exhibit a strongly non-local character It 
would be very interesting to analyze the stationary distri- 
butions obtained here along these lines, although it may 
be impossible to derive exacts analytical expressions as for 
the simple exclusion model. 

Finally, it would be very useful to generalize the Boltzmann- 
Langevin approach of to the computation of the 
complete distribution of the transmitted charge. At this 
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point, we may object that a full quantum-mechanical deriva- 
tion is already available pofl . however, this derivation is 
involving a quenched averaging over the possible impurity 
configurations. For a given mesoscopic sample, this may 
be justified by an ergodic hypothesis, so that averaging 
the transmission matrix over a finite energy window e.V 
becomes equivalent to impurity averaging. But this may 
break down for very small systems. Furthermore, semi- 
classical models are very flexible, since they have been ex- 
tended to treat interaction effects and inelastic processes 
[|ll|[l7||3^,|33|]. Again, we believe a detailed analysis of the 
stationary distribution in configuration space is required 
to make further progress. 

We want to thank O. Verzelen, T. Jolicoeur, G. Bastard 
and R.Ferreira for sharing their computing resources. We 
are also grateful to N.Regnault for some programmer tricks, 
to H. Willaime and P. Tabeling for technical support with 
the preliminary experiment mention in [ ^6| and to H. 
Bouchiat, S. Gueron and B. Reulet for their feed-back. 
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